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ABSTRACT 

We present moderate resolution spectroscopy of extraplanar diffuse ionized 
gas (EDIG) emission in the edge-on spiral galaxy NGC 4302. The spectra were 
obtained with the SparsePak integral field unit (IFU) at the WIYN Observatory. 
The wavelength coverage of the observations includes the [Nil] A 6548, 6583, Ha, 
and [S II] A 6716, 6731 emission lines. The spatial coverage of the IFU includes 
the entirety of the EDIG emission noted in previous imaging studies of this 
galaxy. The spectra are used to construct position-velocity (PV) diagrams at 
several ranges of heights above the midplane. Azimuthal velocities are directly 
extracted from the PV diagrams using the envelope tracing method, and in- 
dicate an extremely steep dropoff in rotational velocity with increasing height, 
with magnitude ~ 30 kms -1 kpc -1 . We find evidence for a radial variation in the 
velocity gradient on the receding side. We have also performed artificial obser- 
vations of galaxy models in an attempt to match the PV diagrams. The results 
of a statistical analysis also favor a gradient of ~ 30 kms -1 kpc -1 . We compare 
these results with an entirely ballistic model of disk-halo flow, and find a strong 
dichotomy between the observed kinematics and those predicted by the model. 
The disagreement is worse than we have found for other galaxies in previous 
studies. 
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The conclusions of this paper are compared to results from two other galaxies, 
NGC 5775 and NGC 891. We find that the vertical gradient in rotation speed, 
per unit EDIG scale height, for all three galaxies is consistent with a constant 
magnitude (within the errors) of approximately 15 — 25kms~ 1 scaleheight -1 , in- 
dependent of radius. This relationship is also true within the galaxy NGC 4302. 
We also discuss how the gradient depends on the distribution and morphology of 
the EDIG and the star formation rates of the galaxies, and consequences for the 
origin of the gas. 

Subject headings: galaxies: halos — galaxies: individual (NGC 4302) — galaxies: 
kinematics and dynamics 



1. Introduction 



Gaseous halos are potentially excellent laboratories for the study of important aspects of 
spiral galaxy evolution. External, and some internal, processes have significant extraplanar 
manif estations. Interactio ns between the galaxy and its surrounding intracluster me dium 
Vollmer et al.ll200ll ). intergalactic medium (such as primordial accretion; e.g 



Oort 



19701 ). and/or neighboring galaxies (e.g.. iMoore et al.lll998[ ). can obviously have effects out- 



side of the disk. Internal processes, too, can ha ye an impact on the environment beyond the 
star forming disk, such as galactic wind s (e.g., Martin 20031: Veilleux et al. 2005 ). and phe- 



nomena described by gal actic fountain ( Shapiro fc Field 197(1 Bregman 198ol ) or chimney 



( iNorman fc Ikeuchilll989l ) models. An understanding of all these processes and their effects 
on galaxy evolution requires an understanding of the properties of the material surrounding 
the disk. 



I n the Milky Way, extraplanar ga s has been discovered both in emission (e.g..lBeuermann et al. 



1985 




Wakker & van Woerden 


Madsen 


2005 


) and absorption 



[e.g., 



Moos et al 



2000 



19981: lLockman fc 



'idoprvhora 



2005 



Savage et al.l 120031 ). It is difficult to 



interpret these observations because our position within the Galaxy makes distance determi- 
nations difficult, and causes confusion along the line of sight; therefore, external systems are 



often targeted to enable study of the entire disk-halo interface. Gas has b een found in t 
los of many external spiral galaxies in the form of neutral hydrogen f e.g.JSwaters et al. 



re ha- 



Matthews fe Wood 



Strickland et al 



199 



0: 



Dettmar 



200 



2003 

T 



Till 



1990l : iRand 



mann et al. 



with extraplana r dust fe.g.. iHowk fe Savage 



Fraternali et al. |2005h . hot X-ray gas fe.g. jBregman fc Pil dis 



20061). and diffuse ionized gas (DIG) (e.g., 



1999 



1996 1: iRossa fc Dettmar 120031: iMiller fc VeilleuxH2003ah . Together 



1997 



1994 



Rand et al. 



Irwin &: Madden! 120061 ) and extraplanar 



star formation (ITullmann et al.l 120031 ) . these observations paint a picture of a multiphase 
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ISM extending well above the star forming disk. 



The origin and evolution of this extraplanar ISM are not yet clear, but two alterna- 
tives aregenerally considered. First, the gas could have been accreted from the IGM (e.g., 
Oortlll970l ). or from companion galaxies (e.g.. Ivan der Hulst fc Sancisil 120051 ) . Alternatively, 
the gas could be participating in a star formation-driven disk-halo flow, such as the one 
described by the fountain model. Lines of evidence which link the morphological and ener- 
getic properties of the extraplanar DIG (EDIG) to the level of star formation activity in the 
underlying disk (see, for example, iRandl Il996l ; iHoopes et al.lll999l ; iRossa fc Dettmarl 120031 ) 
would seem to support the lat ter idea, but clear examples of the former (see the examples in 
van der Hulst Sancisil 120051 ) are observed. It is not yet clear how important the accretion 
process may be in systems with less obvious signs of external influence. To shed light on this 
question, additional evidence may be gained by studying the motion of the halo gas. 

In recent years, several groups have investigated the kinematics of extraplanar gas. The 



H I halos of several ga 



NGC 5775 (Lee et al. 



1997 



Fraternali et al. 



Matthews fc Wood 



axies have been fou nd to rotate more slowl y than the u nderlying disk: 



2001), NGC 2403 flFraternali et alJl200lf ). NGC 891 flSwaters et al. 



20051 ). and possibly the low surface brightness galaxy UGC 7321 
20031) . Early studies of extrapl anar rotation in the optical (e.g., Rand 



2000; iTiillmann et al.ll2000l ; iMiller &: Veilleuxll2003bl ) were limited by the long-slit spectra, in 
the sense that the one-dimensional observations did not allow for disent anglement of the de n- 
sity and velocity information encoded in the line profiles. Recent work ( Heald et al. 2006al lbl. 
hereafter Papers I and II respectively) has benefited greatly from the ready availability of 
integral field units (IFUs), which allow spectra to be obtained simultaneously over a two- 
dimensional portion of the sky. With these observations, the effects of density and velocity 
can be decoupled. In Paper I, TAURUS-II Fabry-Perot observations were used to establish 
that the EDIG in NGC 5775 has a gradient in rotational velocity with height above the 
midplane (z) of ~ 8kms _1 kpc _1 ; in Paper II, analysis of SparsePak observations of NGC 
891 showed a gradient of ~ 15 kms " 1 kpc -1 in the NE quad rant of the halo. The latter result 



agrees with H I results reported by iFraternali et al.1 (120051 ) 



To begin to assess whether the extraplanar gas has an internal origin, the kinematics of 
the EDIG in both N GC 5775 and NGC 891 have been compared to the results of a ballistic 
model (developed by lCollins et al.ll2002l ) of disk-halo flow. Both galaxies are observed to have 
a steeper rot ational velocity gradient th an predicted by the ballistic model (the discrepancy 
described by IFraternali fc Binneyl 12006k for the neutral halo of NGC 891, is in the same 
sense). The difference between observations and model is greater for NGC 891, which has 
a less filamentary EDIG morphology than NGC 5775. In this paper, we address the EDIG 
kinematics of a third edge-on spiral galaxy with well-studied extraplanar gas, NGC 4302. 
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This g alaxy was selected because it is extremely edge-on (see below), and was found by 
Randl (119961 ) to have the most prominent layer of EDIG within his sample of nine galaxies. 
Moreover, the morphology of the EDIG in NGC 4302, which is distinctly different than that 
in either NGC 5775 or NGC 891 (see discussion below and in §ED, enables a comparison of 
EDIG kinematics in morphologically distinct gaseous halos. 

NGC 4302 is classified as Sc in the Third Reference Catalogue of Bright Galaxies (RC3; 
de Vaucouleurs et al.ll 199 ll ). but this classification is somewhat uncert ain because the galax y 
is very nearly edge-on. A bar is likely to be present; see, for example, iLiitticke et al.l (120001 ) . 
NGC 4302 is a member of the Virgo Cl uster, and has a nearby comp anion, NGC 4298, with 
some signs of an interaction (see, e.g., iKoopmann fc Kenneyl 120041 . and discuss i on in § [3]). 
In this paper we assume a distance to NGC 4302 of D = 16.8 Mpc, after iTullyl (119881 ) . To 
estimate the inclinati on angle, we closely examined the Two Micron All Sky Survey (2MASS; 
Skrutskie et al.l 120061 ) J and H images of NGC 4302, both of which show a clearly defined 
dust lane. Based on the apparent offset of the dust lane from the major axis, we estimate 
an inclination angle % ~ 89°, and adopt that value for the remainder of the paper. Small 
deviations from this value will not significantly alter the results of our paper. The dust lane 
is very pronounced in the opt ical bands, and indivi dual features are observed far from the 
midplane (up to ~ 1.5 kpc; see lHowk fc Savagdll999l ). The possible effects of dust extinction 
are discussed where appropriate throughout the paper. A summary of galaxy parameters 
for NGC 4302 is displayed in Table ffl 

The EDIG morphology of NGC 4302 is smoother than that of both NG C 5775 and 



NGC 891. It has been observed as a part of several EDIG imaging surveys (e.g.. iPildis et al. 



1994] ; lRandlll996l ; iRossa fc Dettmarll2000l ). Taken together, the survey images demonstrate 
the existence of an extremely faint, smooth EDIG layer reaching up to 2 ~ 25" (2 kpc), 
and with an apparently sharp radial cutoff at R' ~ 50" (4kpc)§ A single faint plume, 
extending to z ~ —9" (—0.73 kpc; we take z < on the east side of the disk), near 
the nucleus of the galax y was first reported by Pildis et ali (jl994T j. and later confirmed by 
Rossa fc Dettmarl (120001 ) . We are not aware of any other distinct EDIG features (in contrast 
to the extraplanar dust extinction, which has a quite intricate filamentary appearance). We 
note that a n extended dust distrib ution can significantly alter the appearance of the Ha 
emission, as iHowk & Savagd (120001 ) point out in their investigation of NGC 891. However, 



in the case of NGC 4302, we consider it unlikely that the dust extinction has significantly 
altered the appearance of the EDIG in such a way as to result in the observed smooth 



3 To avoid confusion, we use R to represent galactocentric radii, and R' for major axis distance, throughout 
the paper. We express distances in linear units (kpc) in the former case, and offsets in angular units in the 
latter. 
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distribution. 

This paper is arranged as follows. We describe the observations and the data reduction 
steps in §|2J An examination of the halo kinematics is presented in §[31 and the ballistic model 
is compared to the results in §4) In §[5l we compare the EDIG halo kinematics of NGC 4302 
with those of NGC 5775 and NGC 891, and discuss implications regarding the origin of the 
EDIG in these systems. We briefly conclude the paper in §6j 



2. Observations and Data Reduction 



We used the SparsePak fiber array (jBershadv et al l 12004 120051 ) to observe NGC 4302 
during the nig hts of 2004 March 16-17 at the WIYnE 3.5 -m telescope. SparsePak's 82 
4.687"-diameter fibers fed the Bench Spectrograph, which was set up with the 860 lines mm -1 
grating at order 2. This setup yielded a dispersion of 0.462 A pixel -1 and a spectral resolution 
0.60 A (27kms _1 at Ha). The wavelength range covered the [Nil] AA6548, 6583, Ha, 



inst 



and [S II] AA6716, 6731 emission l ines. I n Figure (H the two pointings are shown, overlaid on 
the Ha image of NGC 4302 from iRandl (119961 ) . An observing log is presented in Table [2j 



The po intings were chosen to cover the brightest EDIG emission in the galaxy. The layer 
reported by IRandl (119961 ). defined by \z\ < 25" (2kpc) and \R'\ < 50" (4kpc) was within the 
coverage obtained with our two pointings (see Table [2]). The spectrogra ph setup desc ribed 
above allowed us to detect the EDIG, which is quite faint. As reported by IRandl (119961 ). at a 
height of z = 1 kpc, the EDIG emission measure is 5 — 10 pc cm -6 , about 3 times fainter than 
the EDIG in NGC 891 at the same height. In our spectra, we readily detect the Ha emission 
line at z = 1.3 kpc with a typical peak intensity of 1.76 x 10~ 18 ergs -1 cm" 2 arcsec -2 = 



0.88 pc cm 6 , assuming T = 10 4 K. A higher resolution setup (as was used by iHeald et al. 



2006bl . to observe the brighter EDIG in NGC 891) would have yielded significantly improved 
velocity resolution, but would also have greatly lengthened the exposure times necessary 
to obtain useful signal-to-noise ratios in the velocity profiles at the same heights that we 
consider here. 

The initial data reduction steps were performed in the normal way using the IRAF§ 



4 The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, 
Yale University, and the National Optical Astronomy Observatory. 

5 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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tasks CCDPROC (for bias corrections) and DOHYDRA. This second task performs the flat 
field correction, and calibrates the wavelength scale using observations of a CuAr lamp, 
which were obtained approximately once per hour during the observing run. To estimate 
the accuracy of the wavelength solutions, we calculated the centroids of night-sky emission 
lines using the IRAF task SPLOT. The rms variation in the line centroids was 0.02 A, corre- 
sponding to approximately 0.8 km s _1 at Ha. DOHYDRA also does a rough fiber throughput 
correction using flat field observations, and finally extracts the spectra. At this point, the 
relative aperture throughput has been corrected, but the absolute correction remains un- 
known. Before determining the flux calibration using observations of a spectrophotometric 
standard star (Feige 67), light which fell outside the central fiber because of atmospheric 
blurring of the point spread function or imper fect centering o f the s tar on the fiber must 



be accounted for. To correct for these effects, iBershady et al.l (120051 ) have empirically de- 
termined the appropriate gain corrections (see their Appendix C4) based on measurements 
of the amount of starlight in the central fiber and the six surrounding fibers (which sample 
the wings of the standard star's point spread function). We have measured the amount of 
spillover for each of our spectrophotometric standard observations (obtained approximately 
once per hour throughout the run), applied the appropriate gain correction, and completed 
the flux calibration using the IRAF task CALIBRATE. Based on the uncertainties in the 
gain correction described above, we estimate the accuracy of the flux calibration to be at 
about the 0.7% level; the precision is estimated to be at the 2% level based on the variations 
between the individual observations of Feige 67. 

Once the spectra were fully reduced, the continuum and skyline emission were subtracted 
in the following way. First, short segments of the spectra, containing object emission lines, 
were clipped out. These segments were chosen to be short enough that the continuum 
was well defined and could be removed with a linear fit. Next, sky lines were subtracted 



using our implementation of the "iterative clipping" method described by IBershady et al. 



( 120051 . their Appendix D). This method accounts for the inherent field-dependence of optical 
aberrations in spectrographs. In each wavelength channel, a low-order function (for these 
observations, quadratic) was fit along the slit dimension. In order to obtain a good fit, more 
data points are required than are available by only considering fibers which are entirely free 
of object emission in all wavelength channels. Instead, all 82 fibers are considered together, 
and an iterative clipping scheme is implemented to obtain the correct profile shape along 
the slit dimension (i.e., the slit function for the optical aberrations). This technique is 
best suited for observations where the object emission is Doppler-shifted to widely varying 
wavelengths in the different fibers. In that way, sky lines (which vary little from fiber to 
fiber) may be readily distinguished from object emission and subtracted out. Because of 
the locations of our pointings (one on the receding side and the other on the approaching 
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side), the object emission in our data is not distributed widely in wavelength space. Thus, 
application of this sky subtraction method leads to some oversubtraction of object emission, 
because the iterative clipping routine is unable to reject all of the object signal, leading to an 
overestimation of the normalization of the fit. In the region including Ha and [N II] emission, 
the oversubtraction was minimal; in the region including the [S II] lines, the oversubtraction 
was more severe. We have corrected for this effect by subsequently subtracting the average 
of the (sky-subtracted) spectra corresponding to fibers farthest from the midplane of NGC 
4302 (fibers 6, 17, 19, 36, 38, 39, and 57), which did not receive any object flux. This has 
the effect of correcting the normalization of the fitted slit function, and thereby adding the 
oversubtracted galaxy emission line flux back into the spectra. The technique seems to work 
well for these observations. We tested for the possibility that line shapes were significantly 
affected by the skyline subtraction by comparing Ha, [Nil], and [S II] profile shapes in each 
fiber (because sky lines fall on different parts of the object emission line profiles, the line 
shapes should be affected in different ways for each emission line). This analysis did not 
reveal any systematic line shape errors. Following the sky subtraction, residual continuum 
emission was subtracted from the spectra. 

In § El we will analyze position-velocity (PV) diagrams constructed from the spectra. 
The PV diagrams were constructed by calculating the position, in the galaxy frame, of each 
fiber, and sorting the spectra into the appropriate ranges of R' and z. Each PV diagram 
contains spectra from a different range of z. We choose the separation of pixels in the R' 
axis to match the fiber diameter, thus maximizing the continuity of the diagrams. Note that 
when multiple spectra occupy the same location in a PV diagram, the average is used. 

Because the EDI G in NGC 4302 is faint, with a small scaleheight (~ 550 — 700 pc; 



Collins fc Randl l200ll ). the number of independent PV diagrams that can be generated is 
small. We were able to detect EDIG in PV diagrams created for each of the following ranges 
of z: -20" < z < -5" (—1.6 kpc < z < -0.4 kpc); the major axis -5" < z < 5" 
(-0.4 kpc < z < 0.4 kpc); 5" < z < 20" (0.4 kpc < z < 1.6 kpc); 20" < z < 30 
(1.6 kpc < z < 2.4 kpc). We constructed PV diagrams separately for the Ha line and 
the (equally-weighted) sum of the [NII]A6583 and [S II] A6716 lines. The sum is used to 
increase signal-to-noise in the spectra. The total intensity of [N II] A 6583 + [S II] A 6716 is 
approximately equal to that of the Ha line near the midplane, and increases to almost double 
the Ha in tensity at z ~ 2 kpc as the [Nil] /Ha and [SII]/Ha ratios both increase with z 



v see, e.g. 



Miller fc Veilleuxl l2003bl ) . Note that if there were any line shape errors induced 
by the sky subtraction algorithm, they should appear as systematic differences between the 
results from each type of PV diagram. No such differences are observed. 
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Halo Kinematics 



In this section, we analyze the PV diagrams constructed from the SparsePak spectra 
in an attempt to cha racterize the kinemat i cs of the EDIG in N GC 4302. First, the enve- 
lope tracing method (jSancisi &: Alien! Il979l ; ISofue fc Rubin! 1200 ll ) is used to directly extract 
rotation curves from the PV diagrams. We also present comparisons between the data PV 
diagrams and artificial model galaxy observations. The widths of the velocity profiles along 
the minor axis are examined, and, finally, we investigate radial variations in the observed 
halo kinematics. 



3.1. Envelope Tracing Method 



For an edge-on galaxy such as NGC 4302, care must be take n when deriving rotation 
speeds, for the reasons described by, e.g., iKregel fc van der Kruitl (120041 ) . Briefly, a line of 
sight (LOS) through an edge-on disk crosses many orbits. Each orbit has a different velocity 
projection onto the LOS, and thus contributes differently to the total emission line profile. 
In the case of circular orbits, only one of those projections is maximal: the orbit that is 
intercepted by the LOS at the line of nodes (i.e., the line in the plane of the sky along which 
R = R'). Therefore, only that contribution to the line profile directly provides information 
about the rotation speed. The envelope tracing method uses PV diagrams to extract that 
information and build up rotation curves, under the assumption of circular orbits. To be more 
specific, the rotational (or azimuthal) velocity, V^(R), is found by determining the observed 
radial velocity furthest from systemic in the velocity profile at each R' (after corrections 
for the velocity dispersion of the material, and the velocity resolution of the instrument). 
We consider the results of the envelope tracing method to be valid only at radii where the 
rotation curve is approximately flat. At inner radii, where the derived rotation curve is still 
rising, a rising rotation curve c annot be distinguished from a changing radial density profile 
(see also iFraternali et al.l 120051 ) . 



The details of the envelope tracing algorithm are described in Paper I, and are not 
repeated here. We use the following values of the parameters in eqns. 1 and 2 of Paper 
I. The edge of the velocity profile is located using rj = 0.2 and Ji c = 3cr (where o refers 
to the rms noise), which represent a scaling factor on t he peak intensity of t he line profile 
and a definition of the noise floor, respectively (see also lSofue fc Rubinll200ll ). We also set 
J o| as + of nstr = 40kms _1 (determined using the method described in § I3.2j) . where a gas is 
the velocity dispersion of the DIG. With the last value, we assume that cr gas = 30kms _1 
(recall that o" instr = 27kms _1 ). If this assumption is incorrect, all of the derived rotation 
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speeds will be incorrect by an additive constant, but the derived gradient will be unchanged. 

Because dust extinction in the disk may prevent the velocity profiles from including 
information from the line of nodes, we have also calculated a major axis rotation curve from 
HI observations of NGC 4302 (Chung et al., in prep.). The HI observations were obtained 



configuration; the beam size is 17" x 16" (1.4 kpc x 1.3 kpc) at a position angle —59°, and 
the velocity resolution is <J instT = 10.4 km s -1 . A major axis cut was extracted from the data 
cube and kindly provided to us by J. Kenney; we used the resulting PV diagram to obtain an 
H I rotation curve using the envelope tracing method. An H I velocity dispersion of 10 km s _1 
was assumed. Note that the H I rotation curve was also used to verify the systemic velocity, 
under the assumption that the receding and approaching sides have the same rotation speed 
(~ 175 km s -1 ). The best value of the systemic velocity was confirmed to be 1150 km s^ 1 . 

In Figure [2j we compare the rotation curve derived from the H I data to the one derived 
from the optical emission lines. Note that the radial extent of the H I disk is much greater 
than the SparsePak coverage. The optical rotation curve shows a great deal more scatter, as 
can be expected since the H I beam is much larger than the angular size of a SparsePak fiber. 
Still, the magnitude of the rotation curves appears to be approximately the same, and the 
asymmetry between the shapes of the HI rotation curves on the approaching and receding 
sides seems to be reproduced in the optical. Extinction may explain some of the scatter 
in the optical rotation curve, and the apparently slower rotation speeds derived from the 
optical lines at R > 5 kpc on the approaching (north) side of the disk. There are no other 
obvious signs that extinction may be systematically altering the optical rotation curves. 

Having verified the value of the systemic velocity, and having found no evidence for 
significant extinction errors in the derivation of the rotation speeds, we next utilize PV 
diagrams constructed on the east and west sides of the disk to search for a change in average 
rotation speed with height above the midplane (dV^/dz). In order to maintain consistency 
with previous papers, and to agree with Equation [1] below, we adopt the convention that 
dVrp/dz > for rotation speed decreasing with increasing height. 

For the west side of the galaxy, azimuthal velocity curves were derived for PV diagrams 
constructed from spectra at 5" < z < 20" and 20" < z < 30", and are shown in Figure [3j 
A decrease in rotation speed with increasing height is apparent in the figure, though there 
is a large amount of scatter owing to the faint EDIG emission (the scale height, 7 — 9", 
is considerably less than the positions of the highest fibers, 20" < z < 30" — we are 



with the National Radio Astronomy Observatory (NRAO) 




6 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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observing gas emission up to ~ 4 scale heights from the disk), and perhaps some localized 
extinction. We nevertheless attempt to quantify the decrease in rotation speed with height, 
for the approximately flat part of the rotation curves. By determining the mean azimuthal 
velocity at each height for R > 2kpc (see Figure [3]), and fitting a linear relation to a plot 
of V^-vs.-z, we have obtained values of dV^/dz, which are listed in Table [31 We have also 
tabulated the gradients implied by the change in mean rotation speed from the midplane 
to the range 5" < z < 20" (0.4 kpc < z < 1.6 kpc), and from there to 20" < z < 30" 
(1.6 kpc < z < 2.4 kpc). In order to ensure that the mean azimuthal velocities are not 
unfairly biased by values derived from spectra with low signal-to-noise ratios, we weighted the 
individual azimuthal velocities by the total signal-to-noise in the corresponding spectral line. 
There was no significant difference between the weighted and unweighted means; throughout 
this section, "mean value" refers to the weighted mean. Overall, the magnitude of the 
gradient seems to be approximately 30kms _1 kpc -1 . 

We note that if extinction is preventing us from observing emission from the line of nodes, 
and therefore not truly measuring the azimuthal velocities, then our measured gradients will 
be lower limits. The effects of extinction should diminish with increasing height, meaning 
that the line of sight probes farther into the disk as distance from the midplane increases. In 
a cylindrically rotating halo, such an effect would lead to an apparent rise in rotation speed 
with height. Regardless, comparison with the HI rotation curve indicated that extinction is 
not significantly affecting the derived rotation speeds. 

On the east side of the disk, we are unable to trace the EDIG as far from the plane as 
on the west side. We are able to recover azimuthal velocities in the range —20" < z < —5" 
(—1.6 kpc < z < —0.4 kpc), but no higher. Because of the possibility that a gradient, if 
present, does not begin at z = 0, we hesitate to draw conclusions from the absolute difference 
between azimuthal velocities in the midplane and in the range —20" < z < —5". Instead, 
we compare the azimuthal velocities to the east of the disk to those on the west side, to 
check whether the rotation appears to be consistent on both sides. 

In Figure HI we compare the azimuthal velocity curves. Average azimuthal velocities 
for the approximately flat part of the rotation curves have also been calculated. The plot 
demonstrates very little difference between the rotation of EDIG on the east and west sides 
of the halo at the same absolute distance from the midplane, indicating that the west side 
(the side facing its companion, NGC 4298) does not have significantly different kinematics 
than the opposite side. This does not rule out a connection to a possible tidal interaction, 
but there is no evidence that the gaseous halo on one side of the disk is reacting differently 
than the other. 
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3.2. PV Diagram Modeling 

The envelope tracing method is sensitive to the signal-to-noise ratio of the line profiles 
under consideration, and therefore the apparent gradient found above could be in part due 
to the falling signal-to- noise ratio with increasing height. In order to test for this possibility, 
and to ensure that the density profile of the gas is properly accounted for, we utilize another, 
independent method to derive the kinematics of the EDIG in NGC 4302. This method is 
based on constructing model galaxies with signal-to-noise ratios matched to the data at each 
z, performing artificial SparsePak observations of the models, and comparing modeled PV 
diagrams to the data. 

To generate the galaxy models, we use a modifi ed version of the Groningen Image 



Processing SYstem (GIPSY; Ivan der Hulst et al.l Il992f ) task GALMOD. GALMOD builds 



galaxies from a series of concentric rings. The properties of the rings are specified by: 
the radial density profile of emitting material; the major axis rotation curve; the velocity 
dispersion of the emitting material; the exponential scale height; and the viewing angle of 
the galaxy. Each parameter is allowed to vary from ring to ring, but with the exception of 
the radial density profile we force all parameters to be constant with radius. Note especially 
that by allowing the viewing angle to vary with radius, a warp can be built into the disk. The 
signature of a warp could be confused with a changing rotation speed with height (see, e.g., 



Swaters et al.lll997l ). but we assume that the star forming disk is not warped. Discretization 



noise is included in the model. 

Our modification of the task allows for a linear variation in the rotation speed with 
increasing height, of the form 

v(r a = 1 V ^ R ' Z = ^ for z - Z ^ n 1 

Z) \ V+(R, z = 0) - g [z - z cyl ] for z > z cyl [L) 

where V^(R,z = 0) is the major axis rotation curve, dV/dz is the magnitude of the linear 
gradient (note that dV/dz is positive for a declining rotation speed with increasing height), 
and z cy \ is the height at which the gradient begins. In general, we allow z cy \ to take nonzero 
values, but the present observations do not allow us to set any constraints on the value of 
this parameter, and we choose to leave it set to zero for the remainder of the modeling. The 
vertical (emission mea sure) scale height of the EDIG in NGC 4302 has been measured by 



Collins fc Randl (1200 ll ) to be 0.55 — 0.7kpc; for the models described here, the value 0.65 kpc 
was found to work well. The signal-to-noise ratio was matched to the data by adjusting the 
magnitude of the density distribution at the midplane. 

The output of GALMOD is an artificial data cube. The noise is measured from the 
difference of two models constructed by varying only the random number seed. We have 
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written a script which performs artificial SparsePak observations of the output data cubes, 
by extracting spectra in the pattern that the SparsePak fibers project onto the sky. Once 
the artificial observations have been made, the spectra are handled exactly like the data. 
The size of pixels along the velocity axis in the model is, by construction, the same as in the 
data, and the velocity resolution is set to be the same as in the data. In both the data and 
the model, PV diagrams are constructed separately for the approaching and receding sides 
of the disk and then joined together, except where noted. 

To search for the best model for NGC 4302, we first attempted to estimate the radial 
density profile using the same technique employed in Papers I and II. In those investigations, 
the GIPSY task RADIAL was used to generate radial density profiles by considering intensity 
cuts in Ha images parallel to the major axis. In the case of NGC 4302, models created 
simply using pure exponential disks (scale length 7kpc, the same value used in §HJ), while 
unable to match some of the localized, small-scale features, were overall quite successful in 
matching the data. In retrospect, this result is expected because the EDIG distribution 



is extremely smooth and featureless (see, e.g., iRossa fc Dettmarl 120001 ). We therefore use 
the pure exponential radial density profile to generate our best-fit model. We also use a flat 
rotation curve (175 km s -1 based on the rotation curves derived from the HI data; see Figure 
[2]), and a total velocity dispersion of 40kms _1 (as in the envelope tracing method). 

Next, a sequence of models which vary only in their value of dV/dz was generated. PV 
diagrams were constructed from the models as described above, for the same ranges of z 
used to make the data PV diagrams. Difference PV diagrams are inspected, and the value of 
dV/dz that minimizes the reduced chi-square {xt)i and the value that minimizes the mean 
difference between data and model, are chosen to define the best model for that particular 
data PV diagram. The results are presented in the first six rows of Table HI Although 
we consider the Ha and [N II] A 6583 + [S II] A 6716 PV diagrams separately, in most cases 
the same model is compared with both. The exception is for 20" < z < 30", where the 
increased [Nil] /Ha and [SII]/Ha line ratios require a higher amplitude (by a factor of 1.6) 
for the radial density profile in the [N II] A 6583 + [S II] A 6716 model. Note that for the z- 
range 5" < z < 20" (both in Ha and [Nil] A 6583 + [S II] A 6716) , clumpy structures in the 
data could not be fit well by an axisymmetric model. We therefore mask the three radii 
containing the largest deviations from axisymmetry in the PV diagrams. 

The results of this analysis cluster around a gradient of approximately 30kms _1 kpc™ 1 , 
as did the envelope tracing analysis. There is a hint that the gradient is a bit higher 
than SOkms" 1 kpc -1 . To demonstrate the quality of the PV diagram matching, we present 
overlays in Figures [5] and [7J for the west and east sides, respectively. In the top two rows of 
Figure [5] especially, it is difficult to judge by eye which model provides the best match to the 
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data. Visual comparison is easiest for PV diagrams constructed from spectra at the highest 
z, where the leverage on dV/dz is maximized. The variation in xt with the parameter dV/dz 
is demonstrated in Figure for the models listed in the first six rows of Table HI 

Visual inspection of the overlay plots confirms that a vertical gradient in azimuthal 
velocity of magnitude 30kms _1 kpc _1 provides the best match to the data, but in some 
cases the model with dV/dz = 15kms _1 kpc _1 is also a good match. The conclusion from 
this analysis supports the result of the envelope tracing method, but there is some uncertainty 
in the exact value. 

Because this technique relies on the model density profile being reasonably accurate 
to match the data, we test for the possibility that our adopted exponential disk model is 
biasing the results. To do this, we repeat the procedure described above, but using a flat 
radial density profile. This density profile does not match the data as well as the exponential 
disk, but we were able to analyze the statistics in the same way as before. The results are 
shown in the lower rows of Table HI Despite a large scatter due presumably to the poor 
match between data and model, and an apparent tendency for the best fit gradient to be 
somewhat lower than in the exponential disk models, these results are still consistent with 
dV I dz = 30 km s _1 kpc -1 . 

3.3. Minor Axis Velocity Profile Widths 

Our spectra can also be used to investigate other aspects of the gas kinematics in the 
halo of NGC 4302. Along the minor axis in a galaxy viewed from an edge-on perspective, the 
rotational velocity vectors will all be perpendicular to the line of sight (under the assumption 
of circular rotation). Hence, the only contributions by bulk motions to the width of velocity 
profiles on the minor axis will be made by radial motions. We have examined the spectra 
along (and close to) the minor axis in NGC 4302 to search for signs of radial motions. 
Using the GIPSY task XGAUPROF, Gaussian profiles were fitted to the individual velocity 
profiles (profiles which appeared non-Gaussian were neglected from consideration). The 
velocity dispersions resulting from this fitting procedure were corrected for the instrumental 
resolution (oy = a\ t — of nst ) and are shown in Figure [HJ The plotted errorbars are formal 
errors reported by the fitting routine. There appears to be no evidence for significant change 
in the widths of the velocity profiles as height above the midplane varies. Note that the 
apparent minimum in ay near R' = 0" is expected, because projections of the rotational 
velocity vectors onto the line of sight only vanish on the minor axis. 
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3.4. Radial Variations in dV/dz 

Until this point, we have treated the vertical gradient in rotational velocity as being 
independent of galactocentric radius. Yet an examination of Figure [3] suggests that at least 
on the receding side of the galaxy, the magnitude of t he gradient m ay be variable, changing 



at a radius of about 4kpc. Moreover, as discussed by iRandl (119961 . see especially his Figure 
13), the appearance of the EDIG layer changes at approximately the same radius; the EDIG 
is brighter at small radii (R' < 4kpc), and at outer radii the distribution seems to be more 
dominated by the H II regions. We come back to the relation between halo kinematics and 
morphology in § El for now, we simply seek to justify treating the rotation curves inside and 
outside of R ~ 4 kpc separately. 

To test for a different value of the gradient at different radii on the receding side of 
NGC 4302, we consider separately the rotational velocities displayed in Figure [3] in the 
radial ranges 2.5 kpc < R < 4.25 kpc and 4.25 kpc < R < 6 kpc. (The radial break-point 
was chosen to be the midpoint of the full radial range; an additional test using a slightly 
different break-point, at R = 4.0 kpc, yielded essentially the same results as are presented 
here.) Average rotation speeds at each height were calculated separately in both radial 
ranges. The runs of rotation speed versus height were well fitted with linear relationships, 
and the magnitudes of the fitted gradients are shown in Table [5J Clearly, there is a strong 
difference in the magnitude of the gradient when considering the two radial ranges separately. 
Inside R = 4.25 kpc, the gradient is approximately a factor of 3 shallower than in the outer 
radial range. 

We also performed the same experiment with the rotation speeds on the approaching 
side, which does not appear (from inspection of FigureED to show a significant radial variation 
in velocity gradient. Indeed, the difference between the velocity gradients derived in each 
radial range using the Ha line is within the errorbars. In the case of the rotation speeds 
derived using the sum [N II] A 6583 + [S II] A 6716, there does appear to be a change in gradient 
from the inner disk to the outer disk, but in this case, it is likely to be due to extinction 
lowering the midplane rotation curve in the vicinity of R' ~ 5 kpc (this extinction was also 
invoked in § 13.11 to explain the difference between the H I and EDIG major axis rotation 
curves in Figure [2]). Additionally, separate measurements of the vertical scale height for the 
two radial halves of the approaching side are little different. We therefore conclude that 
there is no evidence for a radial variation in the velocity gradient or the vertical scale height 
on the approaching side, unlike the strong variation observed on the receding side. 
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The Ballistic Model 



The physical explanation for the vertical decrease in azimuthal velocity observed in 
NGC 4302, and other galaxies recently studied, is not yet understood. One possible sce- 
nario is described by the fountain model, which postulates that hot gas is lifted up into 
the halo by star formation activity in the disk. This hot gas would then cool and condense 
into clouds which move through the halo, and eventually rain back down onto the disk. As 
a first step toward understanding the dynamics implied by such a picture, a fully ballistic 



model of disk-halo cycling has b een developed independently by two groups (jCollins et al. 



20021 ; iFraternali fc Binneyll2006l ). Models such as these, which consider the motion of non- 
interacting point masses in a gravitational potential, will be appropriate in cases where the 
density of the cycling clouds is sufficiently greater than that of the surrounding medium, so 
that their motion is essentially unperturbed by hydrodynamics. Whether the relative impor- 
tance of hydrodynamics is indicated by an observable parameter, such as the morphology of 
the EDIG, is unknown, but it seems plausible that filamentary halo structures may imply a 
more ballistic disk-halo flow than smooth diffuse extraplanar gas layers. We return to this 
possibility in § [61 



We have used the ballistic model of 1 Collins et al.l (120021 ) in an attempt to model as 
closely as possible the kinematics of NGC 4302. The model numerically integra tes the orbits 
of ballistic gas clouds in the galactic potential described by lWolfire et al.l (119951 ). The clouds 
are initially placed in an exponential disk, and are launched vertically with an initial velocity 
randomly chosen to be between zero and a maximum "kick velocity," 14 . As the clouds move 
upward and experience a weaker gravitational acceleration, they move radially outward. In 
order to conserve angular momentum, their azimuthal velocity drops as they move outward. 
Thus, a gradient in rotation speed is naturally included in the model. Together with the 
circular speed of the disk (V c , which effectively sets the strength of the galactic potential), 
the kick velocity is the most critical parameter in setting the bulk kinematics of the clouds, as 
well as the scale height of the resulting steady-state halo of clouds. The clouds are assumed 
to have constant temperature, density, and size (and therefore equal Ha intensities); hence 
emission intensity is proportion al to cloud column density along any line of sight. The 
interested reader should refer to ICollins et al.l (120021 ) for a more complete description of the 
model. 

The ballistic model we have considered is only one possible realization of the well known 
"galactic f ounta in" model. The dynamics of this type of circulation was first considered by 
Bregmanl (J1980J) who calculated the trajectory of clouds that condense out of a hot, radia- 
tively cooling halo. The difficulty in constructing such models lies in the lack of knowledge of 
the rotational structure of the gaseous halo that would form from flows of supernova heated 
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gas. Bregman considered two cases, one in which gas flowing into the halo maintains the 
same specific angular momentum as it had in the disk (model "a" in his paper), and another 
in which pressure gradients in the halo lead to cylindrical rotation with the disk (model "b"). 
Our ballistic model corresponds to Bregman's model "a" . The cylindrical model "b" seems 
contraindicated by our observatio nal r esults. More recent models of vertical circulation, e.g. 



de Avillez fc Breitschwerdtl ( 120041 ) and lHarfst et al.l (120061 ) , have a more sophisticated treat 



ment of the interstellar physics, but do not treat a large radial range, nor do they include 
rotation. 

One of the ballistic model inputs is the scale length of the exponential disk. To est imate 



this, w e first inspected a major axis intensity profile extracted from the Ha image of iRand 



( 119961 ). That profile suggested a scale length R sc = 7kpc, and the value was confirmed during 
the analysis described in § 13.21 via comparisons between data and model P V diagrams. The 
circular speed is directly measured from the HI rotation curve in Figure [2j V c = 175 km s -1 . 
The kick velocity is chosen so that the scale height of the model output is clos e to the 
measured scale height of the actual galaxy (550 — 700 pc; ICollins fc Randl 120011 ): Vk = 



60kms _1 . The model is run until the system reaches steady state (after « lGyr). The 
model outputs are the position and velocity of each cloud; these positions and velocities can 
then be used to generate an artificial data cube. 

The azimuthal velocities of the individual gas clouds can be directly extracted from 
the outputs of the ballistic model, and used to generate rotation curves. To compare the 
predictions of the ballistic model to the kinematics of the EDIG in NGC 4302, we have 
extracted rotation curves, using cloud velocities in the same height ranges considered in our 
analysis of the SparsePak data: z < 0.4 kpc, 0.4 kpc < z < 1.6 kpc, and 1.6 kpc < z < 
2.4 kpc. The results are shown in Figure [9j Note that we have extracted rotation curves for 
upward-moving and downward-moving (relative to the disk) clouds separately, in addition 
to considering all clouds together. The first condition would be appropriate if the clouds 
leave the disk as warm, ionized gas, but then cool and return as neutral gas. The second 
condition corresponds to a picture where the clouds leave the disk as hot gas, and return to 
the disk as warm ionized gas. 

Inspection of Figure M reveals immediately that the decrease in rotation velocity with 
height in the ballistic model is extremely shallow. The gradient, calculated in the radial range 
considered during the data analysis described in § 13.11 is approximately 1.1 kms -1 kpc -1 for 
upward-moving clouds, 1.2 kms -1 kpc -1 for downward-moving clouds, and 1.0 kms -1 kpc -1 
for all clouds considered together. These gradients were calculated using only clouds in the 
vertical ranges z < 0.4 kpc and 0.4 kpc < z < 1.6 kpc, because radial redistribution in the 
ballistic model leads to a lack of clouds in the range 1.6 kpc < z < 2.4 kpc for R < 10 kpc, 
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which is a larger radial range than is covered by the SparsePak data. We return to the issue 
of radial redistribution later. 

The gradient in the ballistic model is far less than the value gleaned from the data, 
« 30kms _1 kpc -1 . Recall that the value of the gradient in the ballistic model is driven 
mainly by the value of the ratio of parameters Vk/V c ; in the case of NGC 4302, our best model 
(selected by matching the resulting exponential scale height to the data) is characterized 
by Vk/V c = 60/175 = 0.34. In order to attain a gradient with a higher magnitude, we 
might consider raising the value of the maximum kick velocity. However, a model with a 
kick velocity of 150 km s -1 (V)~/V c = 150/175 = 0.86), while producing a gradient of up to 
« 8kms _1 kpc _1 , also generates a galaxy with a vertical scale height of ~ 31kpc. Clearly, 
increasing the kick velocity will never allow the ballistic model to match the vertical gradient 
in azimuthal velocity measured in the data, while simultaneously matching the scale height 
of the EDIG layer. 

We also note that there is a great deal of flexibility in choosing the shape of the halo 
potential, without drastically affecting the kinematics of the clouds in this model. The reason 
for this, as illustrated in Paper II, is that the region in which the halo potential dominates the 
gravitational acceleration experienced by the clouds is outside of the area typically populated 
by the orbiting clouds. Hence, realistic changes to the halo potential will not dramatically 
change the orbital speeds of the ballistic particles. 

The plots in Figure are missing some data points in the range 1.6 kpc < z < 2.4 kpc 
because of a large amount of radial redistribution predicted by the ballistic model. But is 
the amount of radial redistribution predicted by the model borne out by the data? In Figure 
[TOl we present intensity cuts through the Ha image, as well as through a moment-0 (total 
intensity) map generated from the ballistic model output. The cuts are along the major axis 
and at two heights parallel to the major axis. 

Comparison of the intensity cuts shows that while the distribution of clouds in the 
model at heights 0.4 kpc < z < 1.6 kpc is relatively flat, the intensity distribution in the 
data is much steeper. In the Ha disk, the scale length of the intensity cut is approximately 
7 — 8 kpc; in the range 0.4 kpc < z < 1.6 kpc, the scale length has dropped to about 
2 — 5 kpc, depending on location. Meanwhile, the scale length of the ballistic model intensity 
cuts are ~ 5 kpc in the disk, and about 10 — 20 kpc for 0.4 kpc < z < 1.6 kpc. This shows 
that the outward radial migration in the ballistic model is excessive. The scale length of the 
radial intensity cuts could artificially appear to decrease with height in the data if extinction 
obscures emission from the central regions preferentially at lower z, but we consider this 
explanation unlikely. 
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Discussion 



Now that the ballistic model has been applied to a small sample of spirals with known, 
differing gradients in rotation speed with distance from the midplane, it may be appropriate 
to begin to look for patterns in order to guide future observations. We have summarized 
some important parameters of the galaxies in Table O The listed EDIG scale heights (h z ) 
are electron scale heights, not emission measure scale heights (as are used in § Hj), and were 
determined by fitting (single-component) exponent i al profiles t o Ha images of the three 
galaxies ffrom ICollins et al.l I2OO0I : kand et al.l [l99ol : iRandl Il996l . for NGC 5775, NGC 891, 



and NGC 4302, respectively), over the ranges of R and z indicated in Tabled (the ranges 
are comparable to those used to determine the values of dV/dz, which are also tabulated). 

From the collection of data shown in Table [61 some observed trends are intriguing. 
First, the tendency for the observed velocity gradient to diverge from the ballistic model 
predictions as the EDIG morphology becomes smoother and less filamentary may be an 
indication that smoother halos are intrinsically less governed by ballistic motion. Inversely, 
this lends support to the idea that halos with a more filamentary appearance are better 
described by models of ballistic motion. Effects which have been neglected thus far in our 
ballistic mod el, such as drag, may b e important particularly in the smoother halos. However, 
we note that iFraternali et al.l (120071 ) have recently included drag in their ballistic model (the 
ballistic particles drag against a pre-existing, slowly rotating hot corona); although this 
addition provides a mechanism for removing angular momentum from the ballistic particles, 
thus better matching the HI data for NGC 891, they find that the hot corona "spins up" 
after only ~ 1 Myr, at which point the drag vanishes. 

Another trend which may be suggestive is that the magnitude of the velocity gradient 
is seen to decrease as the filamentary appearance of the halo, the EDIG scale height, and 
the level of star formation activity in the disk (as traced by the surface brightness of FIR 
luminosity, L FIR /-Df 5 ) each increase. This would make sense in a picture where galaxies 
reside in halos populated by a variable mixture of accreting gas and star formation-driven 
fountain-type gas. In such a scenario, the gas associated with accretion, characterized on 
avera ge by low angular moment um, should help to steepen the rotational velocity gradient 
(e.g., IFraternali &: Binneyl 120061 ) ; fountain-type gas, in the absence of such accreting gas, 



results in velocity gradients which are too shallow to reproduce the observed gas kinematics. 
Furthermore, this explanation of the trends we have found would also imply that halos with a 
more filamentary appearance contain proportionally more gas originating in the star forming 
disk, and less accreting low-angular momentum gas. Such halos are associated with disks 
with higher rates of star formation, which is traced here by the quantity Lfir/ 'D% 5 . A larger 
sample of galaxies with varying EDIG morphologies, star formation rates, and observed 
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rotation velocity gradients may be able to refine this picture, and help determine which 
parameters are most closely related. 

Despite the tentative nature of these comparisons, the apparent connection between the 
magnitude of the velocity gradients, and the vertical EDIG scale heights, deserves additional 
comment. The product of the observed value of dV/dz and the EDIG scale height, tabulated 
as dV/dh z , has a large scatter (14 — 36kms _1 scaleheight -1 ), but in each case tends, within 
the errors, toward a value of ~ 20kms -1 scaleheight -1 . If confirmed in a larger sample of 
galaxies, this remarkable relationship should provide an excellent test for theoretical descrip- 
tions of the extraplanar gas (regardless of the details of the physical picture being modeled), 
as it would indicate a kinematic feature common to all gaseous halos. It is comforting that 
the range of EDIG scale heights in our small sample spans roughly a factor of 6, suggesting 
that we have covered a fairly large portion of the parameter space. 

The correlations described above appear to not only hold true for whole galaxies as 
integrated units, but also within individual galaxies. On the receding (south) side of NGC 
4302, the vertical gradient in rotational velocity (31.0 ± 19.8 km s -1 kpc -1 ), expressed in 
units of the scale height (0.56 ± 0.01 kpc), is 17.4 ± 1 1.4 km s -1 scaleheight -1 , consistent 
with the other measurements in Table [6j But on that side of the disk, at a radius of about 
4.25 kpc, the characteristics of the EDIG change (this is clear both from inspection of the 
Ha image displayed in Figured! and from consideration of the measured EDIG vertical scale 
heights). Beyond the same radius, the rotational velocity gradient of the halo gas becomes 
markedly steeper. Therefore, in Table El measurements are tabulated separately for the 
two radial segments of the receding side of NGC 4302 (see also the discussion in § 13. 4p . 
The variation happens in such a way as to keep the velocity gradient, per unit scaleheight, 
approximately constant. 

One explanation for this growing body of obser yational data that seems firmly ruled 



out is the ballistic model. Other types of models (e.g.. iBarnabe et al.ll2006l : iKaufmann et al. 



20061 ) have explained previous analyses of NGC 891, but whether they are able to explain 



the near constan cy of the velocity gradient per unit scaleheight remains to be seen. As 



discussed above, iFraternali fc Binneyl (120061 ) have speculated that the loss of angular mo- 
mentum necessary to match the observations may be provided by an interaction between 
low-angular momentum, accreted gas and the pre-existing halo gas. Whether quantitative 
models of this process will reproduce the observed trends described here also remains to be 
seen. Clearly, more dynamical modeling of extragalactic infall, halo pressure gradients, and 
viscous interactions between the disk and halo would be valuable for interpretation of these 
results. 



A major requirement of any successful model of NGC 4302 is that it not only reproduce 
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the observed constancy of the velocity gradient per unit scaleheight between galaxies, but 
also the constancy of this quantity within a single galaxy with a radially varying scaleheight. 
This provides an intriguing challenge for future interpretation. 



6. Conclusions 

We have presented SparsePak observations of the EDIG emission in the edge-on galaxy 
NGC 4302. By creating PV diagrams from the spectra, using both the Ha line and the 
sum [Nil] A 6583+ [S II] A 6716 independently, we extracted rotation curves at and above the 
midplane with the envelope tracing method, and used those azimuthal velocities to quantify 
the decrease in rotation speed with increasing height. The magnitude of the gradient appears 
to be approximately 30kms _1 kpc" 1 , though there is a good deal of scatter (approximately 
lOkms -1 kpc -1 ). Evidence for a steeper gradient beyond R' ~ 4.25 kpc is found on the 
receding side of NGC 4302. 

As an alternative method for determining the variation in rotation speed with height 
above the disk, we have generated galaxy models constructed using different vertical gradients 
in azimuthal velocity, with signal-to-noise ratios matched to the data at each range of z 
considered. Artificial observations of these models were used to construct PV diagrams, 
which were then compared to the data. Overall, the data tended to favor a gradient consistent 
with the envelope tracing results, dV/dz « 30 km s -1 kpc -1 , but somewhat lower values are 
also possible. 



The ballistic model of I Collins et al.l (120021 ) was used to test the idea that the EDIG in 



NGC 4302 is taking part in a star formation-driven disk-halo flow. In the ballistic model, 
a radial outflow of clouds leads naturally to a decline in rotation speed with height. We 
extracted rotation curves directly from the output of the model, and comparison with the 
data revealed that the predicted gradient is far too shallow. It is possible to obtain higher 
values of the gradient by increasing the kick velocity, but this leads to an extremely large 
vertical scale height in the cloud distribution which is not observed in the data. 

The ballistic model has an additional problem in reproducing the observations. Because 
the clouds undergo large-scale outward radial migrations, the distribution of clouds in strips 
parallel to the major axis rapidly flattens with increasing height, when viewed from an edge- 



on perspective. There is no evidence in the Ha image of iRandl (119961 ) to suggest that this 
radial redistribution is taking place. 



Finally, a comparison was made between the kinematics of the EDIG halos of NGC 4302, 
NGC 5775, and NGC 891. The steepest velocity gradient is measured in NGC 4302. It is 
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intriguing that the kinematics of the halo gas appears to be related to other properties of the 
three galaxies: namely the morphology and vertical scale height of the EDIG halos, and the 
disk-averaged star formation rate (as traced by the quantity Lfir/-D| 5 ). The relationship 
between these parameters is in the sense that would be expected if the halos contain a 
variable mixture of star formation-driven fountain-type gas and low-angular momentum 
accreting gas, and that the steepness of the vertical gradient in rotational velocity is affected 
by the relative proportions of these sources of halo gas. 

It is also very interesting to note that the present observations suggest that the decrease 
in rotational velocity per unit scale height is approximately constant, with magnitude ~ 
15 — 25kms~ 1 scaleheight -1 . This appears to hold true within the receding side of NGC 
4302, where the velocity gradient and vertical scale height both change significantly, but in 
such a way that their product remains approximately constant. 

We thank J. Kenney for providing the H I PV diagram, the observatory staff at NOAO/WIYN 
for their excellent support, and the anonymous referee for a comprehensive report which has 
improved the presentation of this paper. This material is based on work partially supported 
by the National Science Foundation under grant AST 99-86113. 
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Table 1. NGC 4302 Galaxy Parameters 



Parameter 


Value 


R.A. (J2000.0) a 


12 21 42.5 


Decl. (J2000.0) a 


+14 35 54 


Adopted distance b 


16.8 Mpc 


Adopted inclination 


89° 


Position angle 


179° 


Systemic velocity a,c 


1150 kms" 1 


B-band luminosity (L^) d 


4.8 X lO^ergs" 1 


FIR luminosity (-Lfir) 


< 1.3 X lO^crgs" 1 


Star formation rate 


O.SgAfoyr- 1 


Dynamical mass 8 


1.1 X 10 u Af Q 



a From NED. 
b From lTullvl <1988h . 

c Verified using HI data, as discussed in § 13.11 

d Based on the corrected apparent magnitude (m^) given in lde Vaucouleurs et all (1 1 99 lh ■ 

°NGC 4302 and its companion, NGC 4298, are not resolved in the IRAS survey, from which the FIR flux was taken. 
'Derived from LpiR, using the prescription given bv lKennicutj lll998l) . 

s After iHelou et all {l982), but adjusted to reflect the distance (D = 16.8 Mpc) and inclination (i = 89°) adopted in this 
paper. 



Table 2. NGC 4302 SparsePak Observing Log 



Pointing ID 


RA* 


Dccl. a 


Array PA 


R ,h 


z' b 


Exp. Time c 


rms Noisc d 


(see Fig. 1) 


(J2000.0) 


(J2000.0) 


(°) 


(") 


(") 


(hr) 


(erg s _1 cm -2 A -1 ) 


N 


12 21 43.10 


14 36 29.64 


89 


-72 to - 3 


-43 to + 24 


6.3 


3.37 (3.26) X 10~ 18 


S 


12 21 43.10 


14 35 15.96 


89 


+ 2 to + 71 


-42 to + 26 


6.8 


3.22 (3.05) X 10~ 18 



a R.A. and Decl. of fiber 52 (the central non- "sky" fiber in the SparsePak array) for each pointing. 

^Ranges of R f and z' covered by each pointing of the fiber array. "Sky" fibers arc not included in these ranges. R 1 is positive on the south 
(receding) side, and z' is positive to the west. At D — 16.8 Mpc, 12" — 1 kpc. 

c Total exposure time, which is the sum of individual exposures of about 30 minutes each. 

^The rms noise was measured in the continuum near the Ha line for each of the 82 fibers in every pointing. The tabulated values are the mean 
(median) . 
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Table 3. Summary of dV/dz Values for the West Side, using Envelope Tracing Method 



PV Diagram 
(1) 


(2) 


[dVt/dz] ,i 
(3) 


[dV^/dz] h2 
(4) 


App., Ha 

App., [NII]A6583 + [SII]A6716 
Rec, Ha 

Rcc, [Nil] A 6583 + [S II] A 6716 


33 ± 10 
36 ± 8 
29 ± 14 
33 ± 14 


45 ± 27 
35 ± 16 
23 ± 21 
20 ± 21 


24 ± 21 
36 ± 21 
39 ± 28 
61 ± 31 


Note. — (1) PV diagrams were 
Ha line or the sum [N II] A 6583 + 


: constructed from fibers on 
[S II] A 6716. 


the approaching (App.) or receding (Rcc.) side, using either the 



(2) Gradient determined from a linear fit to the average azimuthal velocity at all three heights (— 0.4kpc < z < 0.4 kpc, 
0.4kpc < z < 1.6 kpc, 1.6 kpc < z < 2.4kpc), expressed in units [kms -1 kpc -1 ]. 



(3) Gradient determined from only the average azimuthal velocities at —0.4 kpc < z < 0.4 kpc and 0.4 kpc < z < 1.6 kpc, 
expressed in units [kms - 1 kpc -1 ]. 

(4) Gradient determined from only the average azimuthal velocities at 0.4 kpc < z < 1.6 kpc and 1.6 kpc < z < 2.4 kpc, 
expressed in units [kms -1 kpc -1 ]. 



Table 4. Summary of Determinations of dV/dz using PV Diagram Modeling Method 



PV Diagram 


Emission linc(s) 


Disk model 


[dVj,/dz] x 2 


[dVfi/dzlmean 




(1) 




(2) 


(3) 


(4) 


(5) 


— 1.6 kpc 


< z < 


-0.4 kpc 


Ha a 


Exp 


36 


31 


—1.6 kpc 


< z < 


-0.4 kpc 


NS a 


Exp 


37 


36 


0.4 kpc 


< z < 


1.6 kpc 


Ha 


Exp 


28 


36 


0.4 kpc 


< z < 


1.6 kpc 


NS 


Exp 


31 


31 


1.6 kpc 


< z < 


2.4 kpc 


Ha 


Exp 


27 


22 


1.6 kpc 


< z < 


2.4 kpc 


NS 


Exp 


33 


34 


— 1.6 kpc 


< z < 


-0.4 kpc 


Ha a 


Flat 


23 


34 


— 1.6 kpc 


< z < 


-0.4 kpc 


NS a 


Flat 


23 


15 


0.4 kpc 


< z < 


1.6 kpc 


Ha 


Flat 


22 


38 


0.4 kpc 


< z < 


1.6 kpc 


NS 


Flat 


18 


18 


1.6 kpc 


< z < 


2.4 kpc 


Ha 


Flat 


12 


10 


1.6 kpc 


< z < 


2.4 kpc 


NS 


Flat 


26 


26 



Note. — (1) PV diagrams were constructed from fibers on both the approaching and receding side unless otherwise noted, 
in the listed ranges of z. 

(2) The data PV diagrams were constructed using either the Ha line or the sum [NIIJA6583 + [S II] A 6716 (marked 'NS' here). 

(3) The model PV diagrams were constructed using either an exponential or flat disk, as described in the text. 

(4) Gradient determined by minimizing the xt between data and model PV diagrams, expressed in units [kms -1 kpc -1 ]. 

(5) Gradient determined by minimizing the mean difference between data and model PV diagrams, expressed in units 
[kms -1 kpc -1 ]. 

a In the range —1.6 kpc < z < —0.4 kpc, the spectra on the southern (receding) side were not able to be recreated with a 
realistic density profile, and are therefore not used for this analysis. 
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Table 5. Radial Variation in dV/dz 



PV Diagram 


R.adial range 


[<nVd*] at 


(i) 


(2) 


(3) 


Rec, Ho; 


inner 


22.2 ± 11.6 


Rcc, [NII]A6583 + [SII]A6716 


inner 


23.3 ± 12.9 


Rec, Ha 


outer 


59.5 ± 5.0 


Rcc, [NII]A6583 + [SII]A6716 


outer 


58.8 ± 1.0 


App., Ha 


inner 


33.4 ± 7.2 


App., [NII]A6583 + [SII]A6716 


inner 


43.2 ± 3.1 


App., Ha 


outer 


24.0 ± 3.3 


App., [NII]A6583 + [SII]A6716 


outer 


14.7 ± 3.4 



Note. — (1) PV diagrams were constructed from fibers on the approaching (App.) or receding (Rcc) side, using cither the Hq line or the sum 
[NII]A6583+ [SIIJA6716. 

(2) The radial ranges 2.5 kpc < R < 4.25 kpc (inner) and 4.25 kpc < Ft < 6kpc (outer) were considered separately. 

(3) Gradient determined from a linear fit to the average azimuthal velocity at all three heights ( — 0.4kpc < z < 0.4kpc, 0.4kpc < z < 1.6 kpc, 
1.6 kpc < z < 2.4kpc), expressed in units [kms -1 kpc -1 ]. 



Table 6. Summary of Galaxy Parameters 



Galaxy 


NGC 5775 


NGC 891 




NGC 4302 




Halo Side/Quadrant 


SW side 


NE side 


NE quadrant 


NW quadrant 


SW quadrant 


Observed dV/dz [kms -1 kpc - ^] a 


8 ± 4 


8 ± 4 


17.5 ± 5.9 


34.5 ± 6.4 


22.8 ± 8.7 


59.2 ± 2.5 


i?-rangc [kpc] 


0-12 


0-12 


4-7 


2.5 - 6 


2.5 - 4.25 


4.25 - 6 


z-range [kpc] 


1.2 - 3.6 


1.2 - 3.6 


1.2 - 4.8 


0.4 - 2.4 


0.4 - 2.4 


0.4 - 2.4 


hz [kpc] 


2.2 ± 0.07 


1.7 ± 0.03 


1.40 ± 0.02 


1.05 ± 0.04 


0.63 ± 0.02 


0.35 ± 0.02 


ii-rangc [kpc] 


- 9.6 


- 9.6 


4-7 


2.5 - 6 


2.5 - 4.25 


4.25 - 6 


2-rangc [kpc] 


2.4 - 4.3 


2.4 - 3.6 


0.5 - 2 


0.5 - 1.2 


0.5 - 1.1 


0.5 - 1.1 


dV/dh z [kms 1 scalchcight — 1 ] 


17.6 ± 9.4 


13.6 ± 7.0 


24.5 ± 8.6 


36.2 ± 8.1 


14.4 ± 5.9 


20.7 ± 2.1 


Ball. mod. dV/dz [ km s - 1 kpc - l ] b 


4 


1-2 


1 


iFIR/Das [lO^'ergs -1 kpc -2 ] c 


8.1 


2.2 


< 2.3 d 


EDIG Morphology 


Many bright filaments 


Bright diffuse + filaments 


Faint diffuse 



a From PV diagram modeling in the case of NGC 5775 (Paper I) and envelope tracing analysis in the cases of NGC 891 (Paper II) and NGC 
4302 (this work). 

b Vclocity gradients predicted by the ballistic model for NGC 5775 (Paper I), NGC 891 (Paper II), and NGC 4302 (this work). 
c From lRandl i 19961 ). 

d NGC 4302 and its companion, NGC 4298, arc not well resolved from each other in the IRAS survey, from which the values of Lpjj^ are derived. 
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Fig. 1 . - — (a ) The two pointings of SparsePak, overlaid on the Ha image of NGC 4302 from 



Randl (11996I ). The labels N and S indicate the IDs used for each pointing. The diameter of 
each SparsePak fiber is 4.687 arcsec. (b) The same image, smoothed to 1.5" in order to bring 
out the faint diffuse halo gas. The axes indicate the spatial scale at the adopted distance 
D = 16.8 Mpc. Here and throughout the paper, R! is positive on the receding (south) side, 
and z is positive to the west. 
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Fig. 2. — Comparison of major axis rotation curves derived from the HI (solid line: ap- 
proaching side; dashed line: receding side), Ha (filled squares: approaching side; open 
squares: receding side), and [N II] A 6583 + [S II] A 6716 (filled circles: approaching side; 
open circles: receding side) major axis PV diagrams. Typical errors are approximately 
0.5 — 1 km s _1 for the H I rotation curves (for R < 10 kpc), 0.5 — 2 km s -1 for the Ha rotation 
curves, and 0.5 — 1.5 kms -1 for the [NII]A6583+ [S II] A 6716 rotation curves. 
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Fig. 3. — Plots of azimuthal velocity curves, determined at each height, on the west side of 
the disk. The PV diagrams used to obtain these velocities were constructed from fibers at 
—0.4 kpc < z < 0.4 kpc (filled squares (colored red in electronic edition)), 0.4 kpc < z < 
1.6 kpc (open circles (colored blue in electronic edition)), and 1.6 kpc < z < 2.4 kpc (crosses 
(colored green in electronic edition)). The solid (colored red in electronic edition), dashed 
(colored blue in electronic edition), and dotted (colored green in electronic edition) lines 
respectively indicate the average azimuthal speed determined for each height, as described 
in the text. Azimuthal velocities were derived separately for the approaching side of the disk 
(left panels) and the receding side of the disk (right panels), using the Ha line alone (top 
panels) and the sum [NII]A6583 + [S II] A 6716 (bottom panels). 
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Fig. 4. — Comparison of azimuthal velocities, as derived from PV diagrams constructed 
using spectra from the ranges —1.6 kpc < z < —0.4 kpc (east; filled squares), and 
0.4 kpc < z < 1.6 kpc (west; open circles). Average azimuthal velocities for R > 2.4 kpc 
are plotted as solid and dashed lines for the east and west z-ranges, respectively. Azimuthal 
velocities were derived separately for the approaching side of the disk (left panels) and the 
receding side of the disk (right panels), using the Ha line alone (top panels) and the sum 
[N II] A 6583 + [S II] A 6716 (bottom panels). 
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Fig. 5. — Comparisons between PV diagrams constructed from the data (white contours) 
and the galaxy models described in § 13.21 (black contours). The data PV diagrams were 
constructed from the Ha line (first and third rows) and the sum [Nil] A 6583 + [S II] A 6716 
(second and fourth rows), on the west side, at the height ranges 5" < z < 20" (top pair 
of rows) and 20" < z < 30" (bottom pair of rows). The models were constructed using 
the base parameters described in the text, but different values of the gradient: dV/dz = 
15kms _1 kpc _1 (first column), dV/dz = 30kms _1 kpc _1 (middle column), and dV/dz = 
45kms~ 1 kpc~ 1 (third column). Contour levels for both data and model are 8, 16, 32, 64, 
and 128a in the first and second rows, and 2, 4, 8, and 16<r in the third and fourth rows, where 
a refers to the rms noise. Positive B! is to the south. The systemic velocity is 1150 km s -1 . 
The spatial scale is indicated in the bottom right corner. 
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Fig. 6. — Variation in the xl statistic with the parameter dV/dz, used to select the best 
(exponential-disk) model for PV diagrams constructed from the Ha line (filled squares) 
and the sum [N II] A 6583 + [S II] A 6716 (open circles) at height ranges (a) —1.6 kpc < z < 
—0.4 kpc, (b) 0.4 kpc < z < 1.6 kpc, and (c) 1.6 kpc < z < 2.4 kpc. All xl curves 
have been divided by their minimum value for presentation; the minimum xl values are 
as follows: (a) 4.30 for Ha, 19.1 for [Nil] A 6583 + [S II] A 6716; (b) 57.7 for Ha, 32.2 for 
[Nil] A 6583 + [SII]A 6716; (c) 0.933 for Ha, 2.49 for [N II] A 6583 + [S II] A 6716. The minimum 
xl in each case occurred for the values of dV/ dz shown in Table H] 
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Fig. 7. — Comparison between PV diagrams constructed from the data (white contours) 
and the galaxy models described in § 13.21 (black contours). The data PV diagrams were 
constructed from the sum [Nil] A 6583 + [S II] A 6716, at the height range —20" < z < —5" 
(on the east side of the disk). The models were constructed using the base parameters 
described in the text, but different values of the gradient: dV/dz = 15kms _1 kpc _1 (left), 
dV/dz = 30kms~ 1 kpc~ 1 (middle), and dV/dz = 45kms _1 kpc _1 (right). For the reasons 
described in the text, only data on the approaching side of the galaxy (B! < 0") are consid- 
ered here. Contour levels for both data and model are 4, 8, 16, and 32<r, where a refers to 
the rms noise. North is to the left. The systemic velocity is 1150 km s -1 . The spatial scale 
is indicated in the bottom right corner. 
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Fig. 8. — Velocity dispersions, measured by fitting single Gaussian profiles to spectra 
along and near the minor axis of NGC 4302, and corrected for the instrumental resolu- 
tion (cry = 0"g t — cfnstr)- Dispersions were measured for the Ha line (top panel) and the 
sum [N II] A 6583 + [S II] A 6716 (bottom panel), at three different heights: —5" < z < 5" 
(circles (colored red in electronic edition)); —20" < z < —5" (filled squares (colored blue in 
electronic edition)); 5" < z < 20" (open squares (colored blue in electronic edition)); and 
20" < z < 30" (diamonds (colored green in electronic edition)). 
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Fig. 9. — Plots of azimuthal velocities extracted from the ballistic model. Rotation curves 
were created for clouds in the range z < 0.4 kpc (filled squares), 0.4 kpc < z < 1.6 kpc 
(open circles), and 1.6 kpc < z < 2.4 kpc (crosses). Separate rotation curves were extracted 
for (a) upward-moving clouds, (b) downward-moving clouds, and (c) all clouds. 
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Fig. 10. — Plots of intensity cuts parallel to the major axis in the the Ha image (white lines) 
and the ballistic model (black lines). Profiles are plotted for ranges (a) —20" < z < —5", 
(b) -5" < z < 5", and (c) 5" < z < 20". 



